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ABSTRACT 



The central arcminute of the Perseus cooling flow galaxy, NGC f 275, has been 
mapped with the JCMT in 12 CO(2-l) at 21" resolution, with detections out to at 
least 36" (12 kpc). Within the limits of the resolution and coverage, the distribution 
of gas appears to be roughly E-W, consistent with previous observations of CO, X-ray, 
Ha, and dust emission. The total detected molecular hydrogen mass is ~ 1.6 x 10 10 
Mq, using a Galactic conversion factor. The inner central rotating disk is apparent in 
the data, but the overall distribution is not one of rotation. Rather, the line profiles 
are bluewards asymmetric, consistent with previous observations in HI and [OIII] . We 
suggest that the blueshift may be due to an acquired mean velocity of ~ 150 km s _1 
imparted by the radio jet in the advancing direction. Within the uncertainties of the 
analysis, the available radio energy appears to be sufficient, and the interpretation is 
consistent with that of Bohringer et al. (1993) for displaced X-ray emission. 

We have also made the first observations of 13 CO(2-l) and 12 CO(3-2) emission 
from the central 21" region of NGC 1275 and combined these data with IRAM data 
supplied by Reuter et al. (1993) to form line ratios over equivalent, well-sampled 
regions. An LVG radiative transfer analysis indicates that the line ratios are not well 
reproduced by a single value of kinetic temperature, molecular hydrogen density, and 
abundance per unit velocity gradient. At least two temperatures are suggested by a 
simple two-component LVG model, possibly reflecting a temperature gradient in this 
region. 

Key words: galaxies: individuabNGC 1275- galaxies: clusters: individual: Perseus- 
cooling flows-galaxies:ISM-CO:galaxies- radiative transfer 



1 INTRODUCTION 



Cluster cooling flows present a most intriguing conspiracy. 
X-ray observations have revealed large quantities of hot, dif- 
fuse gas in most if not all clusters (e.g. Fabian 1994). In many 
clusters (70-90%) the central gas densities are high enough 
that the radiative cooling time becomes shorter than the 
Hubble time (e.g. Edge, Stewart, & Fabian 1992; Fabian 
1994). Thus pressure-driven cooling flows are inferred. In 
several clusters, soft X-ray emission has also revealed the 
presence of gas at intermediate temperatures of 10 5 -10 7 K, 
as would be expected for gas condensing out of a cooling 
flow (e.g. Johnstone et al. 1992). If the inferred mass depo- 
sition rates of between 10-1000 Mq are maintained over a 



Hubble time, then ~ 10 - 10 13 M Q of material should be 
deposited in the cores of clusters. 

However, there is little evidence for all of this supposed 
infalling material at other wavelengths. Optical data appar- 
ently rule out star formation with a 'normal' (ie. Galactic) 
mass function (e.g. McNamara & O'Connell 1992). Some 
fraction of the cooling material should exist in cool atomic 
and/or molecular clouds. However, HI emission has not been 
detected in any cooling flow cluster, with upper limits of 
10 9 -10 10 Mq, while HI has been detected in absorption in 
only four clusters (e.g. Jaffe 1992). As well, Annis & Jewitt 
(1993) carried out an unsuccessful search for thermal emis- 
sion from cold dust in 11 clusters, placing upper limits of ~ 
10 10 M Q on the mass of cold gas in these clusters, assum- 
ing a Galactic gas/dust ratio. Searches for molecular gas (ie. 
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CO) have also been unfruitful, with one exception (see be- 
low) in over a dozen clusters (e.g. Grabelski & Ulmer 1990; 
McNamara & Jaffe 1994); again typical upper limits to the 
H 2 mass are ~ 10 10 M . 

The only cooling flow galaxy in which CO has yet been 
detected is the one in the Perseus Cluster (Abell 426) for 
which infall rates of ~ 300-500 Mq per year centered on the 
giant cD galaxy, NGC 1275, have been inferred (Gorenstein 
et al. 1978; Fabian et al. 1981; Mushotzky et al. 1981; Allen 
et al. 1992; Allen & Fabian 1997). Allen & Fabian (1997) 
find excess X-ray absorption from ROSAT data in the cen- 
tral 30" of the cluster implying a total mass of ~ 4 x 10 s 
Mq . A system of Ha filaments extending to w 2' = 40 kpc 
from the core of NGC 1275 has been linked to the cooling gas 
(Lynds 1970; Cowie, Fabian, & Nulsen 1980; Heckman et al. 
1989). NGC 1275 is also home to the highly variable radio 
source 3C84 (e.g. Pedlar et al. 1990). HI, in quantities of ~ 
10 10 Mq, has been detected in absorption against this source 
(Jaffe 1990) and the relativistic outflow appears to be dis- 
placing X-ray emitting gas on arcminute scales (Bohringer 
et al. 1993). A foreground infalling (AVk 3000 km s" 1 ) 
system is also known to exist (Minkowski 1957; Lynds 1970; 
Hu et al. 1983; Caulet et al. 1992) and there is some ev- 
idence that NGC 1275 has experienced a previous merger 
(e.g. Holtzman et al. 1992). Braine et al. (1995) attribute 
the origin of the molecular gas in NGC 1275 to a previous 
merger as do Lester et al. (1995) for the origin of the 100 
/xm emitting dust. 

Previous CO observations of NGC 1275 include single- 
dish detections by Lazareff et al. (1989), and Mirabel, 
Sanders, & Kazes (1989), as well as a partial mapping 
by Reuter et al. (1993) to ~ 20" radius in 12 CO(1-0) 
and 12 CO(2-l) with the Institut de Radioastronomie Mil- 
limetrique (IRAM) telescope. Estimated H2 masses are of 
order Mh 2 ~ 10 10 Mq using Galactic CO to H2 conversion 
factors. High resolution 12 CO (1-0) interferometry has now 
also been carried out by Braine et al. (1995) and Inoue et 
al. (1996). The latter authors find a rotating ring- like CO 
structure with a ring radius of 3.5". A hot (T = 1700 K) 
component of molecular hydrogen has also been directly de- 
tected through 2 /im emission within the central few arcsec- 
onds (Fischer et al. 1987; Kawara & Taniguchi 1993; Inoue 
et al. 1996). No attempt has yet been made to determine the 
physical properties of the molecular gas in this most unusual 
system. 

In this paper, we present 12 CO(2-l), 13 CO(2-l), and 
12 CO(3-2) observations of NGC 1275 using the 15 m James 
Clerk Maxwell Telescope (JCMT). We have mapped the 
galaxy in the 12 CO(2-l) line out to ~ 1 arcmin and in the 
12 CO (3-2) line out to ~ 14 arcsec. In addition, we have ob- 
tained IRAM 12 CO(2-l) and 12 CO(1-0) data, kindly sup- 
plied by Reuter et al. (1993). Together with the JCMT data, 
we perform a line ratio analysis for the central position in 
order to constrain the physical conditions in the molecular 
gas and see whether these are consistent with a cooling flow 
origin. The observations are described in Section 2, the re- 
sults are presented in Section 3, the line ratio analysis is 
discussed in Section 4, some physical quantities are derived 
in Section 5, there is a Discussion in Section 6, and our main 
conclusions are summarized in Section 7. For NGC 1275, we 
take D = 70 Mpc, (cz = 5250 km s -1 ; H = 75 km s" 1 
Mpc -1 ) in which case 1" = 340 pc. 



2 OBSERVATIONS AND DATA REDUCTION 

12 CO(2-l) and 13 CO(2-l) data were obtained during 3 
shifts at the JCMT^] on Mauna Kea during 23-25 Novem- 
ber 1994. The beam-size at 230 GHz ( 12 CO(2-l) ) is 21"and 
is 22"at 220 GHz ( 13 CO(2-l) ), which corresponds to ~ 7 
kpc for Ho = 75 km s _1 Mpc -1 . From scans of Saturn, 
the main-beam efficiency, r) m b, for 12 CO(2-l) and 13 CO(2- 
1) were determined to be 0.69 ± 0.05 and 0.66 ± 0.05 re- 
spectively; these values agree well with that determined by 
JCMT staff between 1 Nov 1994 to Jan 30 1995 (0.69 ± 
0.01; see JCMT WWW page). Unless otherwise indicated, 
all temperatures quoted below are on the main-beam bright- 
ness temperature scale T m (,. The system temperature T sya 
(in T* units) ranged between 500-600K for 12 CO(2-l) and 
600-700K for 13 CO(2-l) . Pointing was done on the central 
radio source 3C84, CRL2688 and Saturn several times per 
night. During an 8 hour run, we found rms drifts in point- 
ing of 3-4 arcsec on average. We chose our (0,0) position to 
be the same as that of Reuter et al. (1993), i.e. a= 03 16 
29.29, 8=+41 19 51.9 (1950), in order to facilitate the line 
ratio analysis. This position is 3. 5" west of the galaxy center 
which is at the 3C84 radio core, i.e. RA(1950) = 3 h 16 m 
29 ?6, DEC(1950) =41° 19' 52" (Pedlar et al. 1990). The 
displacement is << the beam size. 

Integrations were carried out using receiver RxA2 with 
the DAS in its 760 Mhz configuration (~ 1000 km s"\ 0.8 
km s -1 per channel). We beam-switched with a 1 hz chop- 
ping frequency and a 150"throw. CRL 2688 scans were used 
to check the calibration for 12 CO(2-l) and 13 CO(2-l) . Two 
12 CO(2-l) scans of CRL 2688 yield integrated intensities of 
232 K km s -1 (Nov 23, T* units) and 202 K km s -1 (Nov 
25), as compared to the published value of 168 K km s" 1 
on October 20, 1994. For 13 CO(2-l) the measured and pub- 
lished line strengths are 42 and 46 K km s _1 respectively, 
in excellent agreement. Throughout later observing runs, we 
also continued to check that the (0,0) spectrum gave consis- 
tent results. We integrated for 5040 sec at the galaxy center 
in 12 CO(2-l) (Fig. 1), and 5520 sec in 13 CO(2-l) (Fig. 3). 

Subsequent data reduction was done using the UNIX 
SPECX package developed by Rachael Padman. The spec- 
tra were averaged (weighted inversely by T sys ), 20-channel 
binned (giving a final velocity resolution of 16.25 km s _1 ), 
and baseline subtracted (linear baselines were mostly used); 
for 12 CO(2-l) only 3720 sec of data were usable. We also 
carried out a partial mapping in 12 CO (2-1) , with integra- 
tions ranging between 1800-3600 sec for 8 points out to 1 
arcmin from the galaxy center (see below and Section 3.3 
for more details). The sampling was 20" in RA and 15" in 
DEC with beam centers in each row displaced by 10" so that 
full coverage was obtained with the 21" beam, though not 
at the Nyquist rate. 

12 CO(3-2) observations of a 5-point (0,0; -7,+7; 
+7,-7; +7,+7; -7,-7) map at the center of NGC 1275 were 
obtained during service time on July 31 1995. The beamsize 
is 14" at 345 Ghz, and rj m b= 0.525 ± 0.033, as measured 



* The James Clerk Maxwell Telescope is operated by the Royal 
Observatory, Edinburgh on behalf of the Science and Engineer- 
ing Research Council of the United Kingdom, the Netherlands 
Organization for Scientific Research, and the National Research 
Council of Canada. 
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Table 1. Telescope Parameters 



Spectral Line 


Telescope 


Beam FWHM 


Vmb 


12 CO(1-0) 


IRAM 


21" 


0.67 


12 CO(2-l) 


IRAM 


12" 


0.5 


12 CO(2-l) 


JCMT 


21" 


0.69/0.61 1 


13 CO(2-l) 


JCMT 


22" 


0.66 


12 CO(3-2) 


JCMT 


14" 


0.53 



For December 1994 and December 1995 respectively 



from Saturn scans. Receiver RxB3i was used with the DAS 
in its 760 Mhz configuration (~ 600 km s -1 , 0.54 km s _1 
per channel). A total of 3000 sec at 0,0 and 2400 sec at each 
of the other four points was obtained. Pointing done on W3, 
G34.3, CRL618, and Saturn shows ~ 3"drifts over the 8 hour 
shift. Again we have averaged, binned by 20 channels (giving 
a velocity resolution of 11 km s _1 ), and baseline-subtracted. 

Finally, during 7 shifts in December 1995, we extended 
our 12 CO(2-l) mapping by obtaining 12 further points ex- 
tending out to ~ 1 arcmin from the galaxy center; integra- 
tions were typically 5400 sec (see Section 3.3). Receiver A2 
was used in 750 Mhz configuration (0.8 km s _1 per chan- 
nel), with 1 hz beamswitching and a 150"throw. Spectra 
were averaged and 20-channel binned for a final velocity 
resolution of 16.25 km s _1 . Saturn scans were again used 
to determine a r; m 6 value of 0.61 ± 0.05 for 12 CO(2-l) on 
Dec 4/5, where the uncertainty includes known random er- 
rors. The JCMT value is given as 0.69 ± 0.01 (Jan 25/96), 
but the uncertainty in a single measurement is closer to ± 
0.07 (Per Friberg, private comm.) where changes could be 
due to a variety of sources, including sky variations, stand- 
ing waves, uncertainties in beamsize, and errors in pointing 
and focussing (Henry Matthews, private comm.). Thus it is 
important to measure r\ m b during each observing session, as 
we have done. We found the same value of rj mb for each of 
three Saturn scans over Dec 4/5. T sys for 12 CO(2-l) ranged 
between 300-400 K (Dec 3) to 500-600 K (Dec 4/5). Point- 
ing done on W3, 3C84 and Saturn show ~ 3"rms drifts over 
an 8-hour shift. 12 CO(2-l) calibration was checked using 
scans of NGC 1275, which also allowed us to check consis- 
tency with our previous data taken in Dec 1994. Five 600 
sec scans taken at 0,0 yield integrated line strengths agreeing 
with each other to within 25%, and also with our December 
1994 data (e.g. Figure 1) to within 25%. 

The IRAM data were supplied by H.-P. Reuter (private 
communication) and full details can be found in Reuter et 
al. (1993). 

A summary of the spectral lines and telescope parame- 
ters is listed in Table 1. 



3 RESULTS 
3.1 Spectra 

The 12 CO(2-l) and 12 CO(3-2) spectra are shown in Fig. 1 
and 2, respectively. For point (+20,0) of Fig. 1, note that the 
low velocity end of the baseline has been truncated. This is 
because several of the individual scans at this position were 
centered at an offset velocity, necessitating a shift in the 
velocity coordinate before averaging. The non-overlapping 
baseline regions were then truncated, but the integrity of 



the spectral line itself (based on inspection of all individual 
scans with full baselines) has been retained. The 13 CO(2- 
1) spectrum for the center position is shown in Fig. 3 with 
the central 12 CO(2-l) spectrum superimposed. We consider 
the isotopic line to be a marginal detection (cf. the similarly 
placed peaks and troughs in comparison to 12 CO (2-1) ), but 
given the low S/N, we place little weight upon this particular 
spectrum in subsequent analysis. 

The strong radio source, 3C84, occurs at the (0,0) po- 
sition (Pedlar et al. 1990) which has a core strength in the 
range 2 - 6 Jy (Braine et al. 1995; private communication 
to Lester 1995; Inoue et al. 1996) at 115 GHz and a total 
size less than 10 mas (Readhead et al. 1983), or area filling 
factor < 2 x 10~ 7 in our 21" 12 CO(2-l) beam. An inspec- 
tion of the original spectra before and after binning/average 
shows no compelling evidence for absorption features above 
the noise level. This is confirmed by Braine et al. (1995) and 
Inoue et al. (1996), neither of whom detect any 12 CO(1-0) 
absorption in much smaller interferometric beams of 2" and 
4", respectively. 

The integrated line intensities at the measured posi- 
tions, J Tmb dV , are given in Table 2, along with the peak 
Tmb and FWHM. The quoted uncertainties take into ac- 
count the noise in the spectrum and, where possible, uncer- 
tainties in baseline subtraction (note that for most spectra, 
both authors independently computed these values). Cali- 
bration error is not included, since it is less than these other 
uncertainties. The peak Tmb and FWHM were determined 
by fitting Gaussians to the line profiles. For those positions 
where no line was detected, 3cr upper limits are given for 
the peak T m b and J T m b dV, based on the observed scatter 
in the spectra. 

3.2 The 12 CO(2-l) Map and Velocity Field 

Our most extensive observations were obtained in the 
12 CO(2-l) line and a map of integrated intensities is shown 
in Fig. 4a. The map was made by integrating each spec- 
trum over velocity to create a grid of j T mb dV and then 
interpolating onto a finer grid using Langrange polynomials 
using the Grenoble Image and Line Data Analysis System 
(GILDAS) software package. Real detections are observed 
out to 36 arcsec (= 12 kpc) from center. 

The main features in Figure 4a are an east-west elon- 
gation, a peak ~ 15 arcsec east of the galaxy center, and a 
smaller peak near the center itself. The east- west elongation 
has been observed in previous maps, i.e. the single-dish data 
of Reuter et al. (1993) and the interferometer data of Inoue 
et al. (1996). The appropriate comparison of our results is 
with previous single-dish CO(2-l) data, i.e. that obtained 
by Reuter et al. at a resolution of 12" (their Fig. 4). They 
find three peaks, one at (0,0), one roughly coinciding with 
our strongest peak to the east, and another ~ 13" to the 
North- West. After smoothing their data to match our 21" 
beam, the IRAM data shows a peak about 5" to the west of 
centre. Thus, there is some discrepancy between data sets in 
the sense that the JCMT data show a stronger peak to the 
east of center whereas the IRAM data show a stronger peak 
to the west of centre. We cannot reconcile the data sets us- 
ing different baseline subtractions. Possible reasons for the 
discrepancy could include the different beam-size, coverage, 
and sampling between the two data sets. Thus the general 
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Figure 1. 12 CO(2-l) spectra at the various observing positions. The temperature scale is in T m (,. We show only detections and suspected 
detections. 



east-west nature of the emission is confirmed but the relative 
intensities of the east-west peaks requires further observa- 
tions. 

There is good overall agreement between the CO emis- 
sion, the low-velocity Ha emission, the dust, and the hot 
X-ray gas. All share an East-West extension and have a simi- 
lar spatial extent. Like our 12 CO(2-l) distribution (Fig. 4a), 
both the X-ray emission (Bohringer et al. 1993) as well as 
the low- velocity Ha emission (Unger et al. 1990) also display 
peaks on the eastern side of the nucleus. The radio map of 
Pedlar et al. (1990) with ~ 5" resolution, is oriented at 
roughly 90 degrees to the CO and Ha emission. As already 
mentioned, the radio emission seems to have displaced the 
X-ray emitting gas (cf. Figure 1 of Bohringer et al., where 
the ROSAT HRI image is overlaid on the radio contours); 
we return to this point in Section 5.2 and the Discussion. 

Figure 4b shows a crude representation of the veloc- 
ity field in the 12 CO(2-l) line, where we have measured 
the velocity at each of the positions in Figure 1, using the 
SPECX f ind-centroid task. Some of the spectra are low 
S/N, and we have indicated this in Figure 4b by enclos- 
ing the velocity in parathenses [e.g. (5140)]; we have used 

to indicate those positions for which no line was de- 
tected. Figure 4b indicates that, not only is the 12 CO(2-l) 
gas blueshifted at the (0,0) point, but all the gas over the 
mapped area is blueshifted by 100-150 km s~ relative to 
the systemic velocity of the galaxy. (The same conclusion is 
reached if we instead plot the peak velocities using gaussian 



fits to the spectra or whether we try different baseline sub- 
tractions.) This is in contrast to Reuter et al. (1993) who 
found evidence for rotation from their 12 CO (2-1) velocity 
field (their Figure 5), but from their individual 12 CO(2-l) 
spectra (their Figure 2) , it is not clear if this rotation is sig- 
nificant. Our 12 CO(2-l) data show little evidence for rota- 
tion on the larger scales mapped. Since the central 13 CO (2- 
1) spectrum (Fig. 3) shows two peaks, consistent with the 
rotational components observed by Inoue et al. (1996), we 
find that any rotation, if present, appears to be confined to 
a central ring. We discuss this further in Section 6.2. 



4 THE MOLECULAR GAS: LINE RATIOS 
AND EXCITATION CONDITIONS 

4.1 Formation of Line Ratios 

The IRAM 12 CO(2-l) (beam size = 12") and 12 CO(1-0) 
(beam size = 21") J T m bdV quantities of Reuter et al. (1993) 
have been sampled in a 7" grid in the central regions and 
an 11" grid for the outer points. Their calibration error is 
estimated to be « 10% or lower for 12 CO(1-0) and « 10 
- 20% for 12 CO (2-1) (Reuter, private communication). We 
interpolated these data using GILDAS in the same way as 
the JCMT data (see §3) and reproduced the maps of Reuter 
et al. (their Fig. 3 and Fig. 4). The data were then smoothed, 
as required, to a 21 arcsec gaussian beam. 

For the JCMT 12 CO(3-2) data (beam size = 14"), the 
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Figure 2. 12 CO(3-2) spectra immediately around the center 



sampling was at 7 arcscc intervals covering a 28" region. 
These data were also interpolated and smoothed. For these 
data, it was also possible to form a line ratio over a smaller 
14" beam, since the IRAM 12 CO(2-l) 12" resolution map 
could be smoothed to 14". 

Since there was some redundancy in the observations 
(i.e. 12 CO(2-l) for both IRAM and the JCMT), we could 
also check on the reliability of the absolute calibration. To 
do this, we smoothed the interpolated 12" resolution IRAM 
12 CO(2-l) map to 21" and directly compared the emission 
within 21" centered at (0,0) to the JCMT 12 CO(2-l) value 
at (0,0) in a 21" beam. The result is excellent (possibly 
fortuitous) agreement to within 5%. 

The resulting line ratios for the (0,0) position are listed 
in Table 3, where we use the nomenclature, n CO(i-j) = 
f T m b[ n CO(i-j)]dV. The final error bars in the ratio are of 
order 15% to 50% and attempt to take reasonable account of 
noise in the spectra, differences introduced by subtracting 
different baselines or choosing different integration limits, 



differences introduced by interpolation/sampling, and dif- 
ferences (where measured) between telescopes. In the case 
of 13 CO(2-l) , we also include a measured variation between 
the ratio of integrated intensity and ratio of peak intensity 
(Fig. 3). 

4.2 Discussion of Line Ratios 

Our measured 12 CO(2-l) / 12 CO(1-0) ratio of 0.74 ± 0.11 
(Table 3) for the cD galaxy, NGC 1275, falls within the 
range normally observed for spiral galaxies. For example, 
the mean for optically selected nearby spirals is 12 CO(2- 
1) / 12 CO(1-0) = 0.89 ± 0.06, with higher values (~ 1.1) 
measured for perturbed galaxies (Braine & Combes 1992; 
Braine et al. 1993). Values of 0.93 ± 0.22 are found for IR 
bright galaxies (Aalto et al. 1995) and > 1 for some starburst 
galaxies (Loiseau et al. 1990; Braine et al. 1993). Wilson, 
Walker, & Thornley (1997) find values ranging between 0.5 
to 1.07 for specific molecular cloud regions in M33. 
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Figure 3. 13 CO(2-l) spectrum at the (0,0) position (heavy line) with the 12 CO(2— 1) spectrum (light line) superimposed. The apparent 
increase in emission at the low velocity end of the spectrum is due to baseline curvature. 



The isotopic ratio of 12 CO(2-l) / 13 CO(2-l) = 6 ± 3 
(based on our marginal detection of 13 CO (2-1) ) is also typ- 
ical of giant molecular clouds (GMCs), for example, values 
of 4 - 14 with a median of ~ 5 - 6 are found for this ratio 
in the Galaxy (Sanders et al. 1993) and Wilson, Walker, & 
Thornley (1997) find an average of 7 ± 1 for 7 molecular 
clouds in M33. Higher ratios tend to be found for starburst 
galaxies (e.g. 12 for M 82; Tilanus et al. 1991) and IR-bright 
and luminous mergers (> 20 for the CO (1-0) isotopic ratio; 
Aalto et al. 1995). 

In contrast, our 12 CO(3-2) / 12 CO(2-l) of 1.25 ± 0.25 
(21" beam) is higher than normally observed in spirals. 
Typically, the 12 CO(3-2) / 12 CO(2-l) ratio is lower than 
12 CO(2-l) / 12 CO(1-0) . For example, Sanders et al. (1993) 
find 12 CO(2-l) / 12 CO(1-0) and 12 CO(3-2) / 12 CO(2-l) ra- 
tios of ~ 0.8 (ranging from ~ 0.7 to 0.95) and ~ 0.4 (ranging 
from ~ 0.35 to 0.7), respectively, for Galactic GMCs within 
the solar circle. Mean values for the nuclear region of IC 342 
are 1.0 and 0.47 (Irwin & Avery 1992). 

Thus our line ratios for the giant cD galaxy, NGC 1275 
(Table 3) are not untypical of those observed for normal 
spiral galaxies, with the possible exception of the 12 CO(3- 
2) / 12 CO(2-l) ratio, which appears to be somewhat high. 
Similarly high values of 12 CO(3-2) / 12 CO(2-l) (i.e. 1.1 to 
1.3) have also been observed by Wall et al. (1993) in a small 
sample of IR bright galaxies. They require a 2-component 
model to explain their data. We explore this possibility in 
the next section. 



4.3 LVG Analysis and Cloud Properties 

For a line ratio analysis, it is usually assumed that the ob- 
served main beam brightness temperature, T m b„(V), in a 
particular spectral line of frequency, v, for a cloud moving 
at some velocity, V, is just the radiation temperature of the 
clouds, T Rv (V) (see Kutner & Ulich 1981), diluted by a fill- 
ing factor: 

T mK (V) = fUiV) * T Rv (V) [1] 

where the filling factor, ff v (V), includes the area filling fac- 
tor for clouds at velocity, V, and a filling factor in velocity 
(due to dilution within channels) if required. 

In order to make use of the observed (integrated) line 
ratios (Table 3) , it is necessary to assume that the emission 
from each line at velocity, V, is coming from equivalent re- 
gions within the beam, i.e. the filling factors for each line 
at velocity, V, are equivalent. This assumption also ensures 
that the ratios of integrated line strength (which have higher 
S/N) are equal to the ratios of peak T R (V) which is the mod- 
eled quantity. 

We could check this assumption for the two cases 
for which the original (unsmoothed) beam sizes were 
identical. The first is the JCMT 12 CO(2-l) and IRAM 
12 CO(1-0) data. A measurement of the peak T m6 12 CO(2-l) 
/Tmb CO (1-0) ratio gives a result which agrees with the in- 
tegrated ratio (Table 2) to within 5%, well within the quoted 
errors. The second is the JCMT 13 CO(2-l) and 12 CO(2-l) 
data. Since the 13 CO(2-l) line may be more optically thin 
than the 12 CO(2-l) line, this line ratio is more likely than 
the others to suffer from differing filling factors. This in- 
deed appears to be the case, given the relative symmetry of 
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Figure 4. (a) Map of integrated intensities in 12 CO(2— 1) . Contour levels are at 1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 13, and 14 K km s — 1 . 
Structure around the extreme periphery is artificial due to interpolation into regions containing no data, (b) Crude Velocity Field in 
12 CO(2-l) . The velocity has been determined by finding the velocity centroids of each spectrum. Values enclosed in '( )' have more 
uncertainty, and ' — ' indicates that there is no line detected at that position. 



the 13 CO(2-l) line with respect to the systemic velocity, in 
contrast to the observed asymmetry of the 12 CO(2-l) line 
(§3.2). In this case, therefore, the error bar in the 12 CO(2-l) 
/ 13 CO(2-l) line ratio takes into account the variation be- 
tween the ratio of peak T m j, and the ratio of integrated line 
strengths. 

Thus, within the uncertainties quoted and the limita- 



tions of what can be checked, the assumption of equivalent 
filling factors appears to be valid. 

The line ratios of Table 3, applicable to a 21" beam, 
were considered in the context of the Large Velocity Gra- 
dient (LVG) approximation (Goldreich & Kwan 1974; De 
Jong, Chu, & Dalgarno 1975), following Irwin & Avery 
(1992) and using the collisional rates of Flower & Launay 
(1985) which are tabulated from 10 K to 250 K. An 11- 
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Table 2. Line Strengths at each Position 



Table 4. Best-Fit One-Component LVG Results 
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Table 3. CO Line Ratios 
Line Ratio 

12 CO(3-2) / 12 CO(2-l) 
12 CO(3-2) / 12 CO(2-l) 
12 CO(2-l) / 12 CO(1-0) 
12 CO(2-l) / 13 CO(2-l) 



Beam FWHM 
14" 
21" 
21" 
21" 



Value 
1.0 ± 0.15 
1.25 ± 0.25 
0.74 ± 0.11 
6.0 ± 3.0 



level CO molecule is assumed; the excitation temperature of 
the highest level is 305 K. Given a kinetic temperature, Tk, 
molecular hydrogen density, tih 2 , and fractional CO abun- 
dance, X, per unit velocity gradient, (dV/dr), X (dV/dr) -1 , 
as inputs, the radiation temperature, T_r, optical depth, r, 
and excitation temperature, T exc , can be predicted for each 
line, as well as the line ratios. The input quantities can then 
be adjusted until the observed line ratios are reproduced. 

This method makes as few assumptions as possible 
about the physical conditions in the gas (in comparison to, 
say, an LTE model), an approach which is particularly im- 
portant for a cooling flow galaxy in which physical condi- 
tions may be atypical. This approach does, however, assume 
that the molecular clouds within the beam have equivalent 
physical properties, i.e. the beam-averaged line ratios pro- 
vide us with mean beam-averaged physical properties of all 
clouds. We require only an isotopic abundance ratio which 
we take to be I = [ 12 CO ]/[ 13 CO ] = 50 ± 30 (Langer & 



Parameter 

Tfei„ (K) 

riH2 (cm" 3 ) 

X (dV/dr)- 1 (km s" 

I 

T fl [ 12 CO(3-2) ] (K) 



With 13 CO(2-l) 
30 
1 x 10 3 



pc 



1 



Tfl[ 13 CO(2-l) ] 
T«[ 12 CO(1-0) ] 
r[ 12 CO(3-2) ] 
r[ 12 CO(2-l) ] 
r[ 13 CO(2-l) ] 
r[ 12 CO(l-0) ] 
T MC [ 12 CO(3-2) 
T e:cc [ 12 CO(2-l) 
T MC [ 13 CO(2-l) 
T ea; c[ 12 CO(l-0) 



(K) 
(K) 



(K) 
(K) 
(K) 
(K) 



12 CO(3-2) / 12 CO(2-l) 
12 CO(2-l) / 12 CO(1-0) 
12 CO(2-l) / 13 CO(2-l) 



x 

x 10" 

80 
16.5 
20.6 

2.7 
22.7 
42.3 
29.9 

1.5 
10.1 
23.9 
26.0 

8.0 
26.2 
0.80 
0.91 

7.6 



Without 13 CO(2-l) 
20 
3 x 10 3 
1 x 10~ 4 

12.5 
14.6 

16.4 
148 
130 

48.7 
19.7 
19.8 

19.9 
0.86 
0.89 



Penzias 1990; also van Dishoeck & Black 1988). A velocity 
model of the individual clouds is also required. In practise, 
the choice of velocity model is not critical, given the un- 
certainties of the data; here we assume uniformly collapsing 
clouds (cf. Richardson 1985). 

We searched through parameter space over a range of 
density from 10 3 to 10 8 cm~ 3 and abundance per unit ve- 
locity gradient from 10~ 1() to 10~ 2 (km s _1 pc -1 ) -1 , both 
in steps of 10 n ^ 2 , n an integer, kinetic temperature from 10 
K to 250 K in steps of 10 K, and for isotopic abundance 
ratios of 20, 50, and 80. The best-fit (from a \ 2 test) one- 
component results are listed in Table 4, column 2. Since the 
13 CO(2-l) line was only a marginal detection, we then re- 
peated this process, excluding this line. These results are 
listed in Table 4, column 3. 

From Table 4, we find similar results, regardless of the 
inclusion of the isotopic line ratio, except for r. The results 
suggest the presence of cold (20 - 30 K) clouds of low (10 3 
cm~ 3 ) density. A comparison between the predicted line ra- 
tios of Table 4 (last 3 rows) and the observed values of Table 
3, however, shows that even the best-fit result does not re- 
produce the observed 12 CO(3-2) / 12 CO(2-l) and 12 CO(2- 
1) / 12 CO(1-0) ratios to within the quoted errors. Again, if 
the isotope is excluded, the agreement is still poor. We can- 
not exclude the possibility that systematic errors may be 
present, causing us to underestimate the uncertainties of Ta- 
ble 3. However, given the unusual environment of NGC 1275 
(e.g. the presence of the radio core, and the possibilities of 
the cooling flow and past interaction), and considering the 
beam size (21" = 7.14 kpc), we feel it is more likely that the 
clouds in the centre of this galaxy cannot be described by a 
one-component model. 

A simple two-component model could be described in 
which two types of clouds, A and B, are allowed, each with a 
unique value of n(H2), X (dV/dr)~ , Tkin and filling factor. 
Equation [1] could then be replaced by: 



ffA*T RAv + ff B *T RBv 



[4] 



and the line ratios then computed as before. We have indeed 
developed such an LVG model which can accommodate two 
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components. However, this two-component model now has 
7, rather than 3, free parameters and we no longer have 
sufficient data to solve for these parameters. We have, nev- 
ertheless, searched through parameter space (where we have 
allowed the filling factor of component B compared to A to 
range in steps of 10 n ^ 4 , n an integer from -25 to 0) to see if 
there are regimes which can be excluded. 

With a two-component fit, it is possible to match the 
observed line ratios quite well within error bars. As ex- 
pected, the results for the two-component model are gen- 
erally unconstrained and solutions are found over the en- 
tire range of density, temperature and filling factor. How- 
ever, the one conclusion that can be reached from the two- 
component analysis is that it is not possible to find a solution 
for which the two cloud components have the same tempera- 
ture. Varying the density, filling factors, and/or abundance 
per unit velocity gradient between the two components with 
a single temperature cannot reproduce the observed line ra- 
tios. (This conclusion remains, if the isotopic ratio is ex- 
cluded from the analysis.) Our best (but non-unique) two- 
component model which shows the least difference in kinetic 
temperature has a cold (10 K) Component A at a density 
of 1.0 x 10 3 cm" 3 and a hot (170 K) Component B at a 
density of 3.1 x 10 4 cm -3 with the filling factor of the hot 
component about half that of the cold component. Thus the 
simplest interpretation of the results is that temperature 
gradients occur in the molecular gas within the central 21" 
of NGC 1275. This will be discussed further in Section 6. 

5 QUANTITY AND PRESSURE OF THE 
MOLECULAR GAS 

5.1 Total Molecular Gas Mass 

We have estimated the total molecular (H2) gas in the usual 
way, using a Galactic conversion between CO intensity and 
Mh 2 : 

Mh 2 = 5.82 x 10 6 /co(i-o)(7r/4)d 2 , 

where d is the telescope beam diameter in kpc at the 
source distance (here d=7.1 kpc for D=70 Mpc within our 
21" beam). The numerical factor is taken from Sanders, 
Solomon, & Scoville (1984), and is applicable for the 
CO(l-O) transition. We follow McNamara & Jaffe (1994) 
and use the above equation for other transitions by multiply- 
ing the numerical factor by the appropriate line ratios taken 
from Table 3. Our results can be summarized as Mg 2 = 2.8 
x 10 9 M Q ± (8-10) x 10 8 M for the 12 CO(1-0) , 12 CO(2- 
1) , and 12 CO (3-2) transitions, where uncertainties in both 
the distance and the line ratios have been accounted for; no 
uncertainty has been assumed in the Sanders, Solomon, & 
Scoville CO(l— 0) numerical factor. A similar result is found 
for the weaker 13 CO(2-l) transition, but with a correspond- 
ingly larger uncertainty. The total mass over all emission 
mapped (shown by Fig. 4a, Table 2) is Mh 2 = 1.6 x 10 10 
M . 

We can compare our Mh 2 mass of ~ 3 x 10 9 Mq within 
our central 21" beam directly with that of Lazareff et al. 
(1989), and Reuter et al. (1993) who observed NGC 1275 
with the 30m IRAM telescope in the 12 CO (1-0) transition 



with the same beamsize, and found Mh 2 = 6 x 10 9 Mq (all 
comparisons are made assuming a distance of 70 Mpc and 
adjusted to our conversion factor) and 2.7 x 10 9 Mq, respec- 
tively. The latter authors find a total mass of Mh 2 = 8.6 
x 10 9 M© over their total mapped region (a total area of 
~ 40" x 30" which is smaller than that shown in Fig. 5). 
Mirabel, Sanders, & Kazes (1989) obtained 12 CO(1-0) data 
with the 12m NRAO telescope, and found Mh 2 = 3.2 x 10 9 
M© within their 55" beam. Finally, Inoue et al. (1996) find 
Mh 2 of 8.6 x 10 9 and 4.3 x 10 10 M within diameters of 
9"and 65"respectively. Thus, our Mh 2 mass is comparable 
with most previous single dish measurements, though the 
interferometer data give higher results. 

We must stress that the above determinations of the 
molecular gas mass have relied upon a Galactic conversion 
between Mh 2 and CO integrated intensity. Conditions in 
cluster cooling flows are likely to be very different from those 
of Galactic molecular clouds, and the conversion factor may 
depend on such factors as metallicity and the intensity of 
the UV radiation field (e.g. Wilson 1995; Arimoto, Sofue, 
& Tsujimoto 1996). Further observational and theoretical 
work is needed to determine the conversion factor over the 
full range of environments encountered (early and late type 
galaxies, starbursts, cooling flows, etc.). 



5.2 Gas Pressure 

Unfortunately we cannot compute the pressure of the molec- 
ular gas with certainty because we have only rough con- 
straints on the gas temperature and density from our LVG 
analysis (§4.5). However, we can use the values from our 
best two-component solution (Section 4.3) to provide a sam- 
ple solution. Taking T fe4n = 170 K and n« 2 = 3 x 10 4 cm" 3 
for the hot, dense Component B (the pressure from Com- 
ponent A can be ignored), we find a molecular gas pressure 
P ~ 7 x 10" 1() erg cm -3 . Since Component B's density is 
reasonably well constrained and its tabulated temperature 
is a lower limit (§4.3), this result should represent a rough 
lower limit, in the context of the two-component LVG model. 
Bohringer et al. (1993) quote a pressure for the hot X-ray 
emitting gas in the central region of Perseus of 2.5 x 10" 10 
erg cm" 3 , while C. Peres (private communication) finds P 
~ 2.5 x 10" 9 erg cm" 3 at a radius of 8", both from ROSAT 
HRI data (and both values corrected to Ho=75). The lat- 
ter number is the best to compare with our value, and the 
agreement is reasonable given the uncertainties in both our 
LVG analysis and the X-ray determinations near the cluster 
center. We can also compare with the pressure found for the 
optical [SII] emission lines by Heckman et al. (1989): P op t=7 
x 10" 10 erg cm" 3 between 2-5 kpc (6-15"). Thus, there 
is good agreement between the gas pressures as determined 
from X-ray, sub-mm, and optical data. 

Bohringer et al. (1993) overlaid their ROSAT HRI im- 
age of Perseus with the 332 MHz radio data of Pedlar et 
al. (1990), and found that the X-ray gas and radio plasma 
are anti-correlated. In particular, the outer radio lobes are 
coincident with minima in the X-ray surface brightness, and 
Bohringer et al. speculated that the radio plasma has dis- 
placed the thermal gas. Similar results are found for other 
clusters with lobe-dominated radio sources, including Abell 
1795, 2029, 2597 and 4059 (see Sarazin 1997 for references), 



© 1998 RAS, MNRAS 000, 000-000 



10 Terry J. Bridges and Judith A. Irwin 



Abell 2634 (Schindler & Prieto 1997), and Abell 2199 (Owen 
& Eilek 1997). 



6 DISCUSSION 

6.1 Displacement of the Molecular Gas by the 
Radio Plasma? 

We indicated in Section 5.2 that the radio lobes and thermal 
gas have similar pressures, and that the radio lobes and X- 
ray gas avoid each other. Following Bohringer et al. (1993), 
we speculate that the radio lobes are thus displacing the hot 
X-ray gas, the optical line-emitting gas, and the molecular 
gas. If the molecular gas on the side of the galaxy near- 
est to us is being pushed outwards, this could explain the 
blueshift of 100-150 km s" 1 we observe in the 12 CO(2-l) 
and 12 CO (3-2) emission. We now show that this idea can 
also work energetically. 

Roughly speaking, the kinetic energy needed to displace 
the ~ 1.6 x 10 10 Mq (x 1.4 to account for heavier elements) 
of molecular gas at a velocity of 150 km s _1 is ~ 5 x 10 57 
ergs; compare this with the total (minimum) energy in the 
radio halo (Component A) of 3 x 10 57 ergs determined by 
Pedlar et al. (1990). The latter radio value applies to a region 
38" x 15" in size and would actually be higher over a region 
which is comparable to the molecular gas component [Fig. 
4(a)] since more radio emission is present on larger scales. A 
simple linear scaling for area gives a modified radio energy of 
~ 16 x 10 57 ergs, or a factor of ~ 3 higher than the molecular 
gas kinetic energy in the same region. Thus (depending on 
the efficiency of energy transfer and given the uncertainties) 
there appears to be sufficient energy in the radio components 
to account for the velocity displacement of the molecular gas. 

Note that it is unlikely that the blueshifted emission 
represents gas on the far side of the nucleus which is infalling 
as a result of the cooling flow because velocities as high as 
150 km s _1 are not expected for quasistatically contracting 
cooling flow gas. 

6.2 What's Special About NGC 1275? 

Why is NGC 1275/Perseus the only cooling flow cluster with 
detected CO emission? The best upper limits on Mh 2 are 
already a few times 10 9 M Q (e.g. McNamara & Jaffe 1994; 
Brainc & Dupraz 1994; O'Dea et al. 1994), and there are 
at least 14 clusters which would have been detected in CO 
if they contained quantities of molecular gas similar to that 
found in NGC 1275/Perseus. As we discussed in the Intro- 
duction, there have been detections of excess X-ray absorp- 
tion in many clusters (White et al. 1991; Allen & Fabian 
1997). While the physical properties of this absorbing mate- 
rial arc not at all known (see discussion in Sarazin 1997), it 
has been argued that much of it might consist of very cold (~ 
3K) and dense molecular clouds. However, NGC 1275 hap- 
pens to have one of the strongest radio cores of any cooling 
flow central galaxy. This radio core may heat the molecular 
gas to detectable temperatures in NGC 1275/Perseus, while 
the gas remains cold in most other clusters. As we have dis- 
cussed already, the radio plasma is certainly affecting the 
thermal gas on larger scales. 



There are other possibilities. Braine et al. (1995) point 
out that "The FIR and CO emission from NGC 1275 cor- 
respond exactly to what is found in gas-rich spirals. Rather 
than a massive cooling flow, the gas may come from accre- 
tion of one or more gas-rich galaxies." There is evidence for 
a past merger or interaction in NGC 1275 from the shells 
and young globular clusters seen in HST images (note, for 
instance, that NGC 1275 is probably suffering an interac- 
tion with the high- velocity system seen in absorption). In 
this case, the molecular gas would be much warmer and 
hence detectable. As well, the gas could be heated up dur- 
ing the merger event itself. Mergers appear to be rare in 
rich clusters at the present epoch, given their large velocity 
dispersions (e.g. Merritt 1985), and this is consistent with 
the non-detections of CO in other cooling flow clusters. A 
merger of a dust and gas rich spiral could also be the source 
of the thermal emission from dust observed in NGC 1275 
(e.g. Lester et al. 1995). 



7 SUMMARY 

NGC 1275 is a remarkable galaxy with properties not seen 
anywhere else. With its central position in the Perseus cool- 
ing flow, its young blue globular clusters, the evidence for a 
past merger, the foreground infalling system, and its molec- 
ular gas component, it it so far unique amongst all galax- 
ies, not just those of its class. Here we have focussed on 
the molecular gas distribution in NGC 1275, since this is 
the only cooling flow galaxy (so far) which displays a cold 
molecular component. Using the JCMT, we have mapped 
the 12 CO(2-l) distribution of this galaxy at a resolution of 
21", detecting gas out to radii of 36" (12 kpc). Within the 
limitations of our spatial resolution and coverage, we find 
an approximate east-west orientation to the emission, con- 
sistent with that found by previous authors in Ha emission, 
X-ray emission, and the dust component. Evidence for the 
inner molecular rotating disk is also seen, consistent with 
NMA interferometric data. However, on broader scales we 
see no evidence for rotation. 

The CO line profiles are consistently bluewards asym- 
metric, both at the central position in both 12 CO (3-2) and 
12 CO(2-l) and over the entire region mapped in 12 CO(2-l) 
, where the S/N is sufficient to measure this quantity. We in- 
terpret this blueshift as being due to kinetic energy imparted 
by the radio jet/outflow in the advancing direction on the 
near side of the nucleus. With an estimate of total H2 mass 
(assuming Galactic conversion factors) of 1.6 x 10 1() M , 
and velocity offset of 150 km s _1 , we find that ~ 5 x 10 57 
ergs is required to impart this kinetic energy. There appears 
to be sufficient energy in the radio source to account for 
this, given the uncertainties. This result is consistent with 
Bohringer et al. (1993) who found that the X-ray gas has 
been displaced by the radio jet. No corresponding redshifted 
component is seen, possibly because the (receding) radio jet 
is less active. 

Together with additional JCMT 13 CO(2-l) data as well 
as 12 CO(2-l) and 12 CO(1-0) IRAM data from Reuter et al. 
(1993), we have formed three line ratios applicable to a cen- 
tral 21" beam. Two of the line ratios appear to be typical of 
cold Galactic clouds, but the third [ 12 CO(3-2) / 12 CO(2-l) 
] is suggestive of gas under more "active", warmer condi- 
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tions. This is confirmed by an LVG radiative transfer anal- 
ysis which indicates that the observed fine ratios cannot be 
reproduced by a singie component, i.e. a singie set of kinetic 
temperature, molecular hydrogen density, and abundance 
per unit velocity gradient. We have therefore developed a 
simple two-component model to describe the data. In the 
two-component analysis, there are insufficient line ratios to 
fully contrain the set of cloud parameters; however, the re- 
sults indicate that the observed line ratios can be reproduced 
if the two molecular components are at different tempera- 
tures. It may be that a temperature gradient exists in the 
central region possibly due to the radio core. 

Finally, while we note that the characteristics of the 
molecular gas in NGC 1275 have been readily related to the 
nuclear radio activity, these observations do not allow us to 
rule out a possible connection with the cooling flow, itself. 
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